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Abstract 

We report experimental results on the acceleration component probability distribu- 
tion function at R\ = 690 to probabilities of less than 10 -7 . This is an improve- 
ment of more than an order of magnitude over past measurements and allows us 
to conclude that the fourth moment converges and the flatness is approximately 
55. We compare our probability distribution to those predicted by several models 
inspired by non-extensive statistical mechanics. We also look at acceleration compo- 
nent probability distributions conditioned on a velocity component for conditioning 
velocities as high as 3 times the standard deviation and find them to be highly 
non-Gaussian. 



1 Introduction 



Recent developments of new experimental techniques have allowed significant 
progress in the Lagrangian description of turbulence [1,2,3,4,5,6]. Although 
this aspect has been investigated theoretically since the beginning of the 
20th century, the extreme difficulty of recording particle trajectories in high 
Reynolds number turbulence restricted experimental investigations. The avail- 
ability of very high speed imaging sensors gives access to the Lagrangian accel- 
eration in fully developed turbulence. It has been shown that the probability 
density function (PDF) of a Lagrangian acceleration component is strongly 
non-Gaussian, i.e., the PDF exhibits extremely large tails [1,3]. This has been 
also observed in direct numerical simulations of fully developed turbulence [7,8] 

These observations resemble features obtained in the field of non-extensive 
statistical mechanics [9,10,11,12]. This opens two promising perspectives: first 
toward a statistical theory of turbulence and second toward a dynamical mod- 
eling of the motion of fluid particles. The latter comes from the theoretical link 
first established by C. Beck [9] between the Tsallis entropy and a stochastic 



Preprint submitted to Elsevier Science 



2 February 2008 



equation for the Lagrangian acceleration. These developments led to a series 
of models for the PDF of Lagrangian accelerations. The multitude of mod- 
els needs to be tested by an accurate comparison with experimental data, 
especially in the low probability tails of the acceleration PDF. 

Here we present new experimental results that have allowed us to improve 
the estimation of the acceleration PDF. First we describe the modifications 
that have been implemented on the previous experiment, then we show the 
acceleration PDF and compare it to the different classes of models. Finally we 
show and discuss some PDF's of acceleration conditioned on velocity which 
are improvements over results presented previously in Sawford et al. [13]. 



2 The experiment 

The experimental set-up has been described in detail in previous publica- 
tions [1,2,3]. The reader should refer to these for detailed information on the 
experiment. Let us recall briefly the main features. The flow under consider- 
ation is of the Von Karman kind: the water is driven by two counter-rotating 
disks, 20 cm in diameter, 33 cm apart, mounted in a cylindrical container. 
The rotation frequency of the disks is regulated at 3.5 Hz for the experiment 
under consideration here. The integral scale of the flow has been shown to 
be 7 cm and the energy dissipation rate to be 1.14 m 2 s -3 [1]. The transverse 
and axial spatial coordinates of small tracer particles (diameter 46 /im, den- 
sity p = l.OQg/cm 3 ) are recorded by two silicon strip detectors. Each detector 
provides a ID recording of the particle trajectories at a sampling rate up to 
70 kHz. The root-mean-square velocity is 0.43 ms _1 , so that the Taylor based 
Reynolds number is 690. The Kolmogorov time is t v = 0.93 ms, the sam- 
pling frequency is F s = 69.9 kHz so that t v corresponds to 65 samples. We 
study here only one transverse component of the acceleration. Even though 
the velocity components have been shown to be anisotropic, the acceleration 
components are close to isotropic at that Reynolds number [1]. 

Compared to the experiment reported earlier, we are now able to record data 
three time faster (in terms of global efficiency, i.e., including transfer rate from 
the digitizer to the computer, repetition rate, etc.) because of a hardware up- 
grade of the data-acquisition computers. This allowed us to obtain more data 
in a reasonable time. The laser used for illumination has also been upgraded 
to a frequency doubled pulsed YAG laser with mean output power as high as 
35 watts, a typical pulse duration of 200 ns and a pulse repetition rate of up to 
70 kHz. This pulsed laser improves the efficiency of the silicon strip detectors 
due to the short duration of the pulses. 

Another difference lies in the signal processing. As stressed in Voth et al. [1], 
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one has to low-pass filter the trajectories to filter the noise in order to have 
access to the acceleration. In our previous work [1,3], the filtering and differ- 
entiation were achieved by fitting a parabola to the data, in a time window 
whose duration had to be appropriately chosen. Here, this procedure has been 
modified to improve the performance of the data processing algorithm. We 
now use a differentiating and filtering kernel k(r) so that the acceleration is 
obtained by convolution a = k-k x, where * is the convolution product and x 
is the spatial coordinate of the particle. We have chosen a Gaussian kernel 

g( T ) = ^=exp (-^2 J . (1) 

V TTW V w J 



To achieve both a low pass filter and the second order differentiation of tra- 
jectories, one simply has to differentiate twice the Gaussian kernel 

Kr) = -l- 3 (% - l) exp f-^) • (2) 



In practice, however, time is discrete and one has to use a finite time win- 
dow, therefore the kernel is truncated. This requires a modification of the 
normalization. The following kernel was used in the data analysis (for discrete 
times): 



fc(r) = A^-l)exp(-^)+5 if ]rKL (3) 
k(r) — otherwise (4) 



where A and B are normalization constants chosen to ensure that k * 1 = 
and k-k t 2 = 2, i.e., the second derivative of a constant is and the second 
derivative of a parabola is 2. Note that, if L equals one time sample 1/F S , one 
recovers the simple 3 points discrete second derivative (independently of the 
value of w). If L ^> w ^> 1/F S one converges to the continuous kernel g when 
F s goes to infinity. The results obtained from the filtering operation are found 
to be in agreement with those obtained by the parabolic fitting procedure. 
Note that due to the truncation, the actual width of the filter is a slightly 
lower that w. 
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Fig. 1. Lagrangian acceleration PDF at R\ = 690. The kernel parameters are w = 16 
samples, L = 1.5w. The solid line is the stretched exponential fit. Insert: linear plot. 

3 Probability density function of Lagrangian acceleration 



3. 1 Experimental estimation of the PDF 



The results presented here required continuous acquisition of data for two 
weeks. The data set consists of 1.7 10 8 points, which include the new mea- 
surements and the data reported in the previous publications [1,3]. The ex- 
perimental probability density function (PDF) is presented in Fig. 1 as dots. 
The new data allows the estimation of the PDF to reach probabilities more 
than one order of magnitude lower than the results in [1,3], i.e., less than 10 -7 . 
Fig. 1 also shows the stretched exponential fit: 



P{a) = Cexp -a 2 / 1 + 



a(3 






M 


a 





(5) 



with p = 0.513 ± 0.003, 7 = 1.600 ± 0.003, a = 0.563 ± 0.02 and C is a 
normalization constant [C = 0.733). This expression is the same as the one 
proposed earlier [1,3], with slightly different values due to improved accuracy 
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of the fit and perhaps also a difference in the Reynolds number. The initial 
values were p = 0.539, 7 = 1.588, a = 0.508 at R x = 970. This stretched 
exponential fit is shown to follow the experimental data very well. 
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Fig. 2. Relative uncertainty on the Lagrangian acceleration PDF at R\ = 690. The 
rms error is estimated by splitting the data set in 6 subsets. 



In Fig. 2, an estimation of the statistical dispersion of the PDF is shown. It has 
been obtained by splitting the whole data set into 6 subsets. The relative rms 
uncertainty is displayed. One observes a minimum of 3% for small values of the 
acceleration, \a\ < lO^ 2 ) 1 / 2 , for which the number of samples is maximum. 
The uncertainty then increases as \a\ increases. It remains lower than 40% for 
|a| < AOia 2 ) 1 / 2 . 

One consequence of our better estimation of the acceleration PDF is that we 
can conclude that the fourth moment converges. Fig. 3 displays the contribu- 
tion to the fourth moment: a 4 P(a). The first observation is that this function 
decreases as \a\ increases. From our fit, the decay is consistent with a stretched 
exponential: 

P(a) ocexp(-|a| a4 ). (6) 
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Fig. 3. Lagrangian acceleration PDF at R\ = 690. The kernel parameters are w = 16 
samples, L = 1.5w. The solid line is the stretched exponential fit. 

This means that the integral / a 4 P(a) converges. The data gives a flatness of 

(7) 



(a 4 ) 

F = "T"Trr~r = 55 ± 8 . 



{*)'< 



and the stretched exponential fit gives a flatness value equal to F — 55.2. The 
possible influence of the particle response time on this value is expected to be 
low as the characteristic time of the particle is 0.12^. The filtering procedure 
has a limited impact because of the low value of the filter width. Changing the 
width from 0.23t v to 0.31t,, decreases the variance by 10% and the flatness by 
15%. The shape of the PDF is unchanged for a/ (a 2 ) 1 / 2 < 25. Thus we expect 
the flatness to be slightly underestimated but do not expect any qualitative 
change in the shape of the PDF. 



3.2 Comparison with models inspired by non-extensive statistical mechanics 



Among the large class of superstatistics [14], four specific examples have been 
proposed to model the observed acceleration PDF. The first one, proposed 
by C. Beck [9] (and referred to as Beck-% 2 here after), is based on the x 2 - 
distribution. One very interesting feature of this analysis is the possibility of 
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making a connection between a statistical approach and a dynamical approach. 
It allows us to build a stochastic equation for the Lagrangian time dynamics 
of acceleration. In the same spirit, A. M. Reynolds [11] and then C. Beck [12] 
proposed a second model (referred as Reynolds and Beck-Log resp.) based 
on log-normal statistics. The approach of Reynolds is much more ambitious. 
He tries to build a full stochastic equation that models the acceleration. This 
model takes into account not only the small time scales dynamics but also 
the large scales. He takes all numerical parameters from experiments or direct 
numerical simulations reported in the literature. Therefore there are no free 
parameters in his model, once the Reynolds number is chosen. Thus the result 
of the model is not a fit to the experimental data. 



Another more classical approach has been proposed by T. Arimitsu and N. 
Arimitsu [10] (referred as Arimitsu). They use the properties of scale invari- 
ance of the Navier-Stokes equation to build a multifractal description of the 
acceleration PDF. 



These four models are displayed in figures 4 and 5 together with the experi- 
mental data. The first obvious observation is that the ^-distribution model 
(with q = 3/2) does not fit the data. The decay of the tails displays qualita- 
tively incorrect behavior. Indeed, for any value of the parameter q, the model 
predicts a power-law decay which is not observed. For q close to 3/2, which 
gives the best fit, the fourth moment diverges which is not the case for the 
experimental data. The modification of the distribution to log-normal, intro- 
duced later by C. Beck, makes the fit qualitatively much better. In particular 
the flatness, for s 2 = 3.0 takes a finite value F = 60.3. The Reynolds model 
does not agree as well with the data as the others. One reason is that there 
are no free parameters in this model. Second, we low-pass filter the acceler- 
ation at the scale 0.23^. This scale is small but not zero. If the temporal 
position signals from the Reynolds model are low-pass filtered at the same 
scale then the overlap with the experimental data is highly improved [15]. 
The multifractal model by Arimitsu also provides a good agreement with the 
experimental results. Concerning the region of low acceleration region, the 
Arimitsu model captures the correct behavior, contrary to the other models. 
Reynolds' and Beck's log-normal models overestimate the zero acceleration 
probability whereas the Tsallis statistics model underestimates it. 



For all these models the departure from the experimental curve for large ac- 
celerations is best viewed in Fig. 5. Most of the proposed models do not peak 
at the same value of acceleration as the experimental curve. Nevertheless most 
of them are very close to the experiment. 
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Fig. 4. Comparison of the experimental acceleration PDF with the models. Dots: 
data, thin dashed line: Beck-x 2 , q = § [9]. Thick dashed line: Beck-Log, s 2 = 3.0 [12]. 
Solid line: Arimitsu, \i = 0.25 and n = 17.1 [10]. Dot-dashed: Reynolds [11]. Insert: 
linear plot. 




Fig. 5. Contributions to the fourth order moment a 4 P(a). Left: linear plot, right: 
loglog plot. Same symbols as previous figure. 
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Fig. 6. (a) Acceleration PDF's conditional on velocity. From top to bottom, con- 
dition: u = 0, 0.45, 0.89, 1.3, 1.8, 2.2, 2.7 and 3.1 u rms . The curves have been shifted 
for clarity by a factor 2, the top one being unchanged, (b) The conditional mean 
acceleration, (c) The conditional acceleration variance. Dots are the experiment, the 
solid line is a 6th order polynomial fit. 

3.3 Conditional probabilities 

Conditional acceleration statistics are important in many second order stochas- 
tic models [11,13,16,17,18]. A very large number of statistics is required to ob- 
tain acceleration PDF's conditional on large velocities since large velocities are 
rare events. The new large data set allowed us to calculate distributions con- 
ditional on larger velocities than those shown in Sawford et. al. [13]. Figure 6 
shows PDF's of the transverse component of acceleration conditional on the 
transverse component of velocity. The conditional probability distributions are 
highly non-Gaussian. They have almost the same shape as the unconditional 
distribution. The conditional acceleration variance increases with increasing 
velocity which seems to be due mainly to a widening of the tip of the distri- 
bution as can be see in in Figure 6. For homogeneous isotropic turbulence the 
conditional mean acceleration should be zero. Departures from zero reflect the 
anisotropy of our flow although DNS of homogeneous isotropic turbulence has 
also shown slight departures from zero [13]. Sawford et. al. present a scaling 
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argument that indicates that the conditional acceleration component variance 
(Af\ui) ~ uf to leading order. A fit to our data in Fig. 6c. shows good agree- 
ment with this. The dependence of acceleration distributions on the velocity 
violates local homogeneity [19] which is an assumption of Kolmogorov 41 the- 
ory. 



4 Conclusion 

Our new set of data helps to discriminate among different classes of models. 
In particular, the model based on the Tsallis entropy is shown to disagree 
with the experimental observations: the acceleration PDF does not display 
any power law decrease and is well fitted by a stretched exponential decay. 
One consequence is that the fourth moment of the acceleration component is 
finite and we observe a value of the flatness close to 55. The models based on 
log-normal statistics proposed by Beck and Reynolds or multifractal analysis 
(Arimitsu) are indeed close to our experimental observations. Since the only 
difference between Beck's two models is the underlying statistics (Tsallis or 
log- normal), the hypothesis of Tsallis statistics must be discarded in order to 
reproduce the observed behavior of the acceleration PDF. The use of a log- 
normal statistics and the physical link of the (3 parameter in Beck's models 
to the dissipation rate are in agreement with the refined 1962 Kolmogorov- 
Obukhov theory of turbulence which is known to describe well the features of 
intermittency. 

We also show that the joint PDF of acceleration and velocity is not gaus- 
sian. The acceleration variance conditioned on velocity is also not constant. 
These results have important consequences on the structure of the stochastic 
equation required to model dispersion by turbulent flows. 



Acknowledgements 

This work has been supported by NSF grant #9988755. We thank Arthur 
La Porta for valuable discussions and writing substantial parts of the data 
analysis code. 



References 

[1] G. A. Voth, A. La Porta, A. M. Crawford, J. Alexander and E. Bodenschatz, 
"Measurement of particle accelerations in fully developed turbulence" , J. Fluid 



10 



Mech., 469 pp. 121-160 (2002). 

[2] G. A. Voth, A. La Porta, A. M. Crawford, E. Bodenschatz C. Ward and J. 
Alexander, "A silicon strip detector system for high resolution particle tracking 
in turbulence", Rev. Sci. Instrum., 72(12), pp. 4348-4353 (2001). 

[3] A. La Porta, G. A. Voth, A. M. Crawford, J. Alexander and E. Bodenschatz, 
"Fluid particle accelerations in fully developed turbulence", Nature, 409 pp. 
1017-1019 (2001). 

[4] N. Mordant, P. Metz, O. Michel and J.-F. Pinton, "Measurement of Lagrangian 
velocity in fully developed turbulence", Physical Review Letters 87(21) 214501 
(2001). 

[5] N. Mordant, J. Delour, E. Leveque, A. Arneodo and J.-F. Pinton, "Long time 
correlations in Lagrangian dynamics: a key to intermittency in turbulence", 
Physical Review Letters 89 (25) 254502 (2002). 

[6] S. Ott and J. Mann, "An experimental investigation of the relative diffusion 
of particle pairs in three-dimensional turbulent flow", J. Fluid Mech. 422 pp. 
207-223 (2000). 

[7] P. Vedula & P. K. Yeung "Similarity scaling of acceleration and pressure 
statistics in numerical simulations of isotropic turbulence", Phys. Fluids, 11, 
1208 (1999). 

[8] T. Gotoh and D. Fukayama, " Pressure spectrum in homogeneous turbulence", 
Phys. Rev. Lett, 86, 3775 (2001). 

[9] C. Beck, "Dynamical Foundations of Nonextensive Statistical Mechanics", 
Phys. Rev. Lett. 87, 180601 (2001). 

[10] T. Aritmitsu and N. Arimitsu, "Multifractal analysis of fluid particle 
accelerations in turbulence", arXiv:cond-mat/0210274 (2002). 

[11] A. M. Reynolds, "On the application of nonextensive statistics to Lagrangian 
turbulence", Phys. Fluids, 15(1) LI (2003). 

[12] C. Beck, "Lagrangian acceleration statistics in turbulent flows", 



[13] B. L. Sawford, P. K. Yeung, M. S. Borgas, P. Vedula, A. La Porta, A. M. 
Crawford, E. Bodenschatz, 'Acceleration Variance and Conditional Variance 
Statistics in Turbulence', to be submitted. 

[14] C. Beck and E. G. D. Cohen, "Superstatistics" , |arXiv:cond-mat / 0205097 
(2002). 

[15] A. M. Reynolds, private communication. 

[16] B. L. Sawford, "Reynolds number effects in Lagrangian stochastic models of 
turbulent dispersion", Phy. Fluids, A3, p. 1577 (1991). 



11 



[17] B. L. Sawford, P. K. Yeung, "Eulerian acceleration statistics as a discriminator 
between Lagrangian stochastic models in uniform shear flow", Phys. Fluids, 12, 
pp. 2033 (2000). 

[18] M. S. Borgas, P. K. Yeung, "Conditional fluid-particle accelerations in 
turbulence", Theoret. Comput. Fluid Dyn. bf 11 pp. 69-93 (1998). 

[19] S. B. Pope, "Turbulent Flows", Cambridge, NY: Cambridge University Press 
2000 p 190. 



12 



